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Abstract 

This note proves that arbitrary local gates together with any entangling bipartite gate V are univer- 
sal. Previously this was known only when access to both V and V' was given, or when approximate 
universality was demanded. 

A common situation in quantum computing is that we can apply only a limited set S C Ud of unitary 
gates to some <i-dimensional system. The first question we want to ask in this situation is whether gates 
from S can (approximately) generate any gate in VUd = Ud/Ui (the set of all d x d unitary matrices up 
to an overall phase). When this is possible, we say that S is (approximately) universal. See [1,3,4,7] for 
original work on this subject, or Sect 4.5 of [9] or Chapter 8 of [8] for reviews. 

Formally, S is universal (for VUd) if, for all W G VUd, there exists U\, . . . ,U k G S such that 

w = u k u k _ x ---u 2 u u 

whereas U is approximately universal (for VUd) if, for all W G VUd and all e > 0, there exists U\ , . . . , U k G 
S such that 

d(W,U k U k . 1 ---U 2 U 1 ) <e. (1) 

Here d(-, •) can be any metric, but for concreteness we will take it to be the VUd analogue of operator 
distance: 

Similar definitions could also be made for Ud, other groups, or even semigroups. 

A natural way to understand universality is in terms of the group generated by S, which we denote (S), 
and define to be smallest subgroup of VUd that contains S. An alternate and more constructive definition is 
that (S) consists of all products of a finite number of elements of S or their inverses. When 5 contains its 
own inverses (i.e. S = S~ l := {x : x" 1 G S}) then (S) provides a concise way to understand universality: 
S is universal iff (S) = VUd and S is approximately universal iff (S) is dense in VUd. 

But what if S does not contain its own inverses? The equivalence between approximate universality and 
(5) being dense in VUd still holds. One direction remains trivial: if S is approximately universal then (S) 
is dense in VUd- The easiest way to prove the converse is with simultaneous Diophantine approximation, 
which implies that for any U G VUd and for any e > 0, there exists n > such that d(U n , U" 1 ) < e. The 
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proof is due to Dirichlet, and for completeness we include it here . For any W G VUd and e > 0, the fact 
that {S} is dense in VUd means that there exists an | -approximation to W of the form U^ 1 . . . Uj^ 1 , with 
each Ui G S. Now we replace each U^ 1 term with U™ 1 for m satisfying ||£7" 4 — U^ l \\ < e/2k. By the 
triangle inequality this yields an e-approximation to W out of a finite sequence of unitaries from S. 

The case of exact universality is more difficult, and is the subject of the current note. Again if S is 
universal then (S) = VUd, and again we would like to argue that the converse holds. Unfortunately this 
statement is not known to be true, and there may well be counter-examples along the lines of the Banach- 
Tarski paradox. However in the special case where S contains a non-trivial one-parameter subgroup then we 
can prove that universality with inverses implies universality without inverses. In fact we prove something a 
little stronger: not only can any element of VUd be written as a finite product of elements from S, but there 
is a uniform upper bound on the length of these products. If we define S L to be the set of products of L 
elements from S, then we can prove 

Theorem 1. 

(a) Suppose S C VUd, (S) = VUd and there exists a Hermitian matrix H such that H is not proportional 
to the identity and e lHt G S for all t G R. Then S is exactly universal for VUd- In fact there exists an 
integer L such that S L = VUd- 

(b) Suppose S C Ud, (S) = Ud and there exists a Hermitian matrix H such that H has nonzero trace, H 
is not proportional to the identity and e lHt G S for all t G R. Then S is exactly universal for Ud, and 
there exists L such that S L = Ud- 

The main interest of this theorem is in its application to the setting of a bipartite quantum system where 
local unitaries are free and nonlocal operations are restricted. Say that d = d^ds and that S = Ud A x 
Ud B U {V}, where U dA x Ud B is embedded in Ud A d B according to (Ua, Ub) — > Ua^Ub and V is some 
arbitrary unitary in Ud A d B ■ In other words, we can perform V as well as arbitrary local unitaries, meaning 
unitaries of the form Ua <S> Ub- Say that V is imprimitive if there exists \^pa) G C dA , \<pb) G C rfs such 
that V(|<^a) ® \pb)) is entangled. Equivalently V is imprimitive if it cannot be written as Ua <8> Ub for any 
U A G U dA , U B G U dB , nor, if d A = d B , as SWAP • (U A ® U B ). Then [1] proved that {S) = VU d if and 
only V is imprimitive. It was claimed in [1] that in fact S was exactly universal when V is imprimitive, but 
their proof assumed that G S. Theorem 1 then fills in the missing step in the proof of [1], and together 
with the fact that local unitaries contain at least one nontrivial one -parameter subgroup and the results of 
[1], we obtain 

Corollary 2. If S = Ud A x U dB U {V} and V is imprimitive then S is exactly universal for Ud A d B - In 
fact, there exists an integer L such that S = Ud A d B - 

This corollary is used in [6] to prove that unitary gates have the same communication capacities with or 
without the requirement that clean protocols be used. Exact universality there is used to show that a protocol 
(possibly inefficient) exists for exact communication using a fixed bipartite unitary gates supplemented by 
arbitrary local operations. Now we turn to the proof of Theorem 1. 

'We prove the claim for U 6 Wj, and the VUd result will follow from the fact that ignoring a global phase can only decrease 
distance. Let the eigenvalues of U be ( e 2 ^ iai , . . . , e 2 * iaid ) for some a £ (R/Z) d . Here (K/Z) d is the d-dimensional torus, which 
can be obtained by gluing together opposite faces of the hypercube [0, l] d . Note that under the Loo-norm, a ball of radius e/2 
will have volume e d . Thus, if n > l/e d then the set {0, a, 2a, . . . , (n — l)a} will have two distinct points, rua and 112a, with 
lima - 11201100 < e. If n' = \n 2 - m| then < n! < n and ||n'a||oo < e. This implies that \\U n _1 - i7 _1 ||oo < |1 — e le | = 
2 sin e/2 < e. 
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Proof. We start with an overview of the proof, and then discuss the details of each step. Let G denote the 
group we are working with, which could be either VUd or Ud, and let m = d 2 — 1 if G = VUd or m = d 2 
if G = Ud- Note that G is an m-dimensional manifold [5]. 

(1) We will define a smooth (i.e. infinitely differentiable) map / from IR m to G. It will have the property 
that dfo (its derivative at the point 0) is non-singular. 

(2) We will construct a map / : W 71 — > G such that dfo is non-singular and there exists an integer I such 
that f{x) € S e for all x e R m . 

(3) We will construct an open neighborhood N of the identity matrix I £ G such that N C S i+i ' for 
some integer £'. 

(4) We will show that G = N n for some integer n, and thus that G = S n ^ +e '\ 



Step 1: For some U\ , . . . , U m € G to be determined later, we define 

f(x) = Uie iHxi ulu 2 e iHx2 ul ■ ■ ■ U m e iHxm Ul 
The partial derivatives at x = are given by 



in • 



dXj 



We would like to choose Ui, . . . , U m so that the UjHUj are linearly independent. Consider first the G = 
VlAd case. Then the space of Hermitian traceless matrices (which we call su^) is a d 2 — 1-dimensional irrep 
of G, so the span of {UHU^ : U € G} is equal to all of su^- Thus, there exists a basis of m = d 2 — 1 
matrices of the form UjHJjj. 

When G = Ud, the tangent space is instead the set of Hermitian matrices Ud, which decomposes into 
irreps as Ud = siy © MI. Since H is neither traceless nor proportional to /, it has nonzero overlap with both 
irreps. Again we would like to show that the span of {UHW : U € G} (which we denote by h) is equal 
to Ud- First, we use the fact that Ud acts transitively on matrices of fixed spectrum. Averaging over all d\ 
diagonal matrices isospectral to H we obtain find that (tr H)I/d, which we have assumed is nonzero, is in 
h. Second, we replace H with H — (tr H)I/d (which is in rj and Slid) and use the result for VUd to show 
that the span of su^ C h. Thus rj equals all of u^. Since t) was spanned by matrices of the form UHU^ , this 
means we can choose a set of d 2 linearly independent matrices U\HUl, . . . , U m HUm to form a basis for 

In either case, dfo has m linearly independent columns of length m, and thus is non-singular. Denote 
the smallest singular value of dfo by e. 

Step 2: Since (S) = G, S is approximately universal and so we can approximate Uj and llj with products 

of elements of S, which we call Uj and Uj respectively. Demand that each approximation be accurate to 
within e/4m, so that the total error is < e/2. We then define / as follows: 



f{x) := U x e inxi UlU 2 e lUx *Ul ■ ■ ■ U m e lHXm U, 



in ■ 



Note that d(/(x), fix)) < e/2 for all x £ M. m . Since e is the smallest singular value of dfo, then the smallest 
singular value of dfo must be > e/2, and thus dfo is non-singular. 
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Additionally, each e lHxj G S and each Uj and U- is a product of a finite number of elements from S, so 
there exists I such that f(x) G 5^ for all x G M m . 

Step 3: According to the inverse function theorem (see e.g. [5]), / is a local diffeomoiphism at 0. This means 
that there exists a neighborhood X of such that f(X) is a neighborhood of /(0) and / : X — > f(X) is 
a diffeomorphism (one-to-one, onto, smooth and such that f^ 1 is also smooth). Let Bs(U) := {V : 
d{U, V) < 5} denote the open ball of radius 5 around U. Since f(X) is a neighborhood of /(0), there 
exists 5 > such that f?2<5(/(0)) C f(X). Now we again use the approximate universality of S to construct 
a 5-approximation to /(0)" 1 , which we call V. Then V ■ f(X) contains B$(I) =: N. Additionally, if 
V G S e ' then AT C V ■ f{X) C S e+e ' . 

Step 4: If n > vr/2 sin- 1 (5/2) then = G. This is because G = {e iH : \\H\loo < vr} (optionally 

modulo overall phase) and B S (I) = {e iH : < 2si n - 1 (V2)}. Thus G = S n ( e+e '\ □ 

We conclude with some open questions. First, it would be nice to know the exact conditions on 5 for 
which (S) = G implies exact universality. A perhaps more important question is that of efficiency. If S 
is approximately universal and contains its own inverses, then the Solovay-Kitaev fheorem[2, 8] states that 
any gate can approximated to an accuracy e by S for £ = poly log(l/e). But if 5 does not contain its own 
inverses, the best bound known on I is the trivial poly(l/e) bound from Dirichlet's theorem. 
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